Fast Relaxation Solvers for Hyperbolic-Elliptic Phase Transition Problems

Phase transition problems in compressible media can be modelled by mixed hyperbolicelliptic systems of conservation laws. Within this approach phase boundaries are understood as shock waves that satisfy additional constraints, sometimes called kinetic relations. In recent years several tracking-type algorithms have been suggested for numerical approximation. Typically a core piece of these algorithms is the usage of exact Riemann solvers incorporating the kinetic relation at the location of phase boundaries. However, exact Riemann solvers are computationally expensive or even not available. In this paper we present a class of approximate Riemann solvers for hyperbolic-elliptic models that relies on a generalized relaxation procedure. It preserves in particular the kinetic relation for phase boundaries exactly and gives for isolated phase transitions the correct solutions. In combination with a novel sub-iteration procedure the approximate Riemann solvers are used in the tracking algorithms. The efficiency of the approach is validated on a barotropic system with linear kinetic relation where exact Riemann solvers are available. For a nonlinear kinetic relation and a thermoelastic system we use the new method to gain information on the Riemann problem. Up to our knowledge an exact solution for arbitrary Riemann data is currently not available in these cases. Preprint Series Issue No. 2011 Stuttgart Research Centre for Simulation Technology (SRC SimTech) SimTech – Cluster of Excellence Pfaffenwaldring 7a 70569 Stuttgart [email protected] www.simtech.uni-stuttgart.de Fast Relaxation Solvers for Hyperbolic-Elliptic Phase Transition Problems Ch. Chalons∗ F. Coquel† P. Engel‡ Ch. Rohde‡